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Charge and spin Josephson currents in a ballistic superconductor-ferromagnet-superconductor 
junction with spin-triplet pairing symmetry are studied using the quasiclassical Eilenberger equation. 
The gap vector of superconductors has an arbitrary relative angle with respect to magnetization of 
the ferromagnetic layer. We clarify the effects of the thickness of ferromagnetic layer and magnitude 
of the magnetization on the Josephson charge and spin currents. We find that 0-7r transition 
can occur except for the case that the exchange field and d-vector are in nearly perpendicular 
configuration. We also show how spin current flows due to misorientation between the exchange 
field and d-vector. 
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I. INTRODUCTION 

Since both of superconductivity and ferromagnetism 
are antagonistic ordered phases of matter, for a few 
decades after discovery of superconductivity, the in- 
terplay between superconductivity and ferromagnetism 
had not been a subject of intensive research interest. 
However, recently the study of interplay between su- 
perconductor and ferromagnet has been revived and, 
in particular, proximity effect and Andreev reflection 
in superconductor-ferromagnet junctions has attracted 
much attentionP"^ 

When a singlet superconductor has placed in proxim- 
ity to a ferromagnetic layer, one can find triplet pairing 
correlations in the ferromagnetic layer and this compo- 
nent of pairing provides interesting effectsP^ On the 
other hand, spin-triplet superconductor have been ex- 
tensively investigated due to its anomalous featuresp2H32] 
Proximity of spin-triplet superconductor and ferromag- 
net is more interesting than the singlet one in the sense 
that both of spin-triplet superconductor and ferromagnet 
have a magnetic nature. 

In RefsPSHSU the authors demonstrated presence of 
both charge and spin currents in the systems consist- 
ing of a thin barrier of ferromagnet sandwiched between 
two triplet superconductors. They obtained a spin cur- 
rent due to the coupling of the ferromagnetic moment 
to the spin of triplet Cooper pairs. Also, they found a 
0-7T transition in the triplet superconductor-ferromagnet- 
superconductor (SFS) Josephson junction generated by 
the misalignment of the d-vectors of the triplet super- 
conductors with respect to the magnetic moment of the 
ferromagnet. In contrast to RefsP*"^ where the ferro- 
magnet is modeled as S function potential, in this paper 
we allow for arbitrary length of the ferromagnet, which 
is a more realistic modeling, based on the quasiclassical 
Eilenberger equation. 

In this paper, we investigate a triplet SFS struc- 
ture with arbitrary misalignment between magnetization 
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FIG. 1: (Color online) Scheme of triplet SFS junction. Here, 
L,R label the left and right half-spaces. Two superconduct- 
ing bulks are separated by a ferromagnetic layer with the 
exchange field h. 



of the ferromagnet and gap vector of superconductors 
(FigfTj), by changing the exchange field and thickness of 
ferromagnetic layer. We find that charge and spin cur- 
rents show 0-7r transition except for the case that the two 
vectors are in nearly perpendicular configuration 

The organization of this paper is as follows. In Secjll] 
the quasiclassical equations for Green functions are pre- 
sented. In Sec |IH[ we calculate Green functions of 
the system analytically, presenting the formulas for the 
Green functions to calculate the charge and spin current 
densities. 



In Sec IV we show the results of spin and 
charge currents. The paper will be finished with the con- 
clusions in SecfVl 
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FIG. 2: (Color online)The x component of the charge cur- 
rent versus the phase difference <f> for I = £/2, T = 0.05Tc, 
and different h and a. Currents are given in units of jo,e = 
feiV(0)7j F A(0). 
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FIG. 3: (Color online) The x component of the charge current 
versus the phase difference cf> for h — 5A, T — 0.05Tc, and 
different / and a. 



II. FORMALISM AND BASIC EQUATIONS 

We consider a clean SFS structure with a homogenous 
ferromagnet of thickness I and bulk triplet superconduc- 
tors(see FigfTj). The ferromagnetic layer is characterized 
by an exchange field h. There is a misorientation angle a 
between the exchange field of ferromagnet and order pa- 
rameters of superconductors (d- vectors) . Interfaces are 
fully transparent and magnetically inactive. The thick- 
ness I is larger than the Fermi wave length and smaller 
than the elastic mean free path. Then, we can use a qua- 
siclassical description in the clean limit, and apply the 



Eilenberger equation in this limilPSSAl 



v^Vg 



£m<5- 3 + i(A + h),<7 



= 0, 



(1) 



with the normalization condition gg = I. Here, e m = 
nT(2m + 1) are discrete Matsubara energies, T is the 
temperature, vp is the Fermi velocity and 03 = T3 (?) I in 
which T3 is the third Pauli matrix in particle-hole space. 
&j (J = 1, 2, 3) denote Pauli matrices in spin space in the 
following. The Matsubara propagator 
in the standard fornr: 



can be written 



!J 



gi + gi • a (g 2 + g2 ■ o") «0" 2 

i&2 (ff3 + g3 • O") 54 - 02 (g4 ' &)<72 



(2) 



where & = (&i, 02, 173) is the vector of Pauli matrices in 
spin space. Also, we can write h as follows: 



h & 







(h.«T) 



(3) 



We use Eilenberger equation for both ferromagnetic 
and superconducting materials. For superconductors we 
set h = 0. The matrix structure of the off-diagonal self 
energy A in the Nambu space is 



A 



(d ■ (r)i&2 

i<72(d* • <t) 



(4) 



In the ferromagnet we set A = and h = 
(0, /isina, hcosa). In this paper, the unitary state, 
d x d* = 0, is investigated. 

Solutions of Eq. (fTl) have to satisfy the conditions for 
Green functions in the bulks of the superconductors: 



9 (±00) 



di.2 = d Ljfl (v>) exp ± 



£m<73 + zAi ;2 

di, 2 | 2 



(5) 
(6) 



where <f> is the external phase difference between the or- 
der parameters of the bulks of superconductors. Eq. ([!]) 
have to be supplemented by the continuity conditions at 
the interfaces between ferromagnet and superconductors. 
For all quasiparticle trajectories, the Green functions sat- 
isfy the boundary conditions both in the right and left 
bulks as well as at the interfaces. 



III. ANALYTICAL RESULTS OF GREEN 
FUNCTIONS 




FIG. 4: (Color online) The normal component (x component) 
of the spin current at the left interface, versus the phase dif- 
ference 0, for I = f/2, T = 0.05T C , different h and a = it/ 4. 
Currents are calculated in units of jo,s = ^hN (0)vfA(0) . 
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step-like function of gap vector which is zero in the fer- 
romagnet and finite value in the superconducting part. 
It has been shown that the absolute value of a self- 
consistently detemined order parameter decreases near 
the interface between superconductor and normal metal, 
while its dependence on the direction of Fermi velocity re- 
mains unaltered^. In RefsPSJEH a qualitative agreement 
between self-consistent and non-self-consistent calcula- 
tions for unconventional Josephson junctions is found. 
Also, it is observed that the results of non-self-consistent 
calculation in RefpS are coincident with the experimen- 
tal data in Ref.33. Also, the results by non-self-consistent 
treatment in RefF^ are similar to the experimental re- 
sults in Ref. 49 . Consequently, we believe that non-self- 
consistent approximation doesn't change results qualita- 
tively. Therefore, in our calculations, a simple model 
of the constant order parameter up to the interface is 
considered to obtain qualitative results. Using such a 
model, the analytical expressions for the charge and spin 
currents in the ferromagnet and superconductors can be 
obtained for a specified form of the order parameter. 
The expre ssion s for the charge and spin currents are the 
followingsPMl 



je(r) 



2iTreTN(0)J2(vFgi(v F ,r,£m)), (7) 



j s< (r) = nhTN (0) V (v F (e, • gl (v F , r, e m ))) (8) 
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FIG. 5: (Color online) The x component of the spin current 
at the left interface, versus the phase difference <f> for h = 5A, 
T = 0.0570, a = tt/4 and different I. 



In principle, Eilenbcrger equations (TTj) should be sup- 
plemented with a self-consistent equation for gap vec- 
tor. While these numerical self-consistent calculations 
can be used to investigate the spatial variation of gap 
vector in the superconductors and also the dependence 
of gap vector on temperature, in this paper we do not 
apply them. In our analytical calculations, we consider a 



where e.;= (x,y,z) and N(0) is the density of states in 
the normal state at the Fermi energy. We consider the 
order parameter as follows: 



d n (T,v F ,r) = rf„(T,v F ,r)z. 



(9) 



Solving Eilenberger equation in ferromagnet region and 
superconductors under the continuity of solutions across 
the interfaces, x — 0,1 and the boundary conditions at 
the bulks, we obtain the Green's functions in the ferro- 
magnet region. The Green's functions are given by 



9lF 

and 



rj{p 2 - 1) (p 2 + 2p[fi + sin 2 a{2/3 2 - /3 - 1)] + l) 



D 
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D 



sina[Bcos(^f) - Asin(f^)] + Ccosa 



where 



A 
B 
C 
D 

P 

V 



n. 



sin 2a[/3 - 1] [p 2 (20 + 1) + 2p(/3 sin 2 a • 



cos 2 a) 



1]. 



sin 2a sia(2/ii/v a? ) [p 2 (2,0- 1) + 2p(/3 sin 2 a + cos 2 a) + 1], 
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(e m +r)Qi)(s m + r/Sl 2 ) 

, + |di, 2 | 2 . 



cos 2 a) 




hi A 



04 






03 






02 






0.1 


' - 






2 4 6 8 


1 






0.4 


\\ 




a = 7i/4 


U 


0.3 











0.2 

0.1 




^c"" v 









1 


2 


3 4 5 



0.4 


\ 


\ 


a= 


= 7l/2 


0.3 






\. 




0.2 






s. 




0.1 














12 3 4 


5 



ft/A 



A/A 



//!; 



11% 



FIG. 6: (Color online) The x component of the critical charge 
current versus h for T — 0.05Tc and different I and different 



FIG. 7: (Color online) The x component of the critical charge 
current versus the thickness of the ferromagnetic layer at T = 
0.05Tc for different h and different a 



For parallel alignment of gap vector and exchange held, 
a = 0, Green's functions reduce to 



9if 

glF 



r](p 2 - 1) 
p 2 + 2pf3 + 1 ' 
2rjp sm(2hl/v x ) . 

p 2 + 2pf3 + 1 



(12) 



It is seen from this expression that Green's functions and 
hence charge and spin currents depend on the exchange 
held through 0. Also, for the perpendicular configuration 



a 



we obtain 



9if = 



r)(p-l) 



1 



glF = 



(13) 



where the Green's function are independent of the ex- 
change field. As a result, we do not see any — 7r tran- 
sition in this case^M 



IV. NUMERICAL RESULTS OF CHARGE AND 
SPIN CURRENTS 

Here, we have considered p— wave superconductors 
with the model of order parameter as 



d(T,v F ,r) = A(T)(k x +ik v ) 



(14) 



which would be realized in S^RuOj 1 ^. The function 
A (T) describes the dependence of the order parameter d 
on the temperature T. We have used Green functions in 



Eqs.( 10 1 and ( 11 ) and calculated charge and spin currents 



for our specified model of the order parameter. For repre- 
sentative values of h, I and a, we have shown charge and 
spin currents in Figs. |2][5]as a function of phase difference 
between superconductors. In the figures, £(T) = J^f T) 
denotes a superconductor coherence length at the bulk of 
superconductor. By changing the parameters, we find a 
tendency towards the 0-7T transition. In particular, we see 
a stronger tendency towards the 0-7r transition at a = 0. 
In our formalism, there are in general three compo- 
nents (s x , s y and s z ) of spin current due to the presence 











FIG. 8: (Color online) The x component of the critical charge 
and spin currents versus a for T = 0.05Tc- Left panels are 
for I = £ and right panels are for I = 2£. 



of the ferromagnetic layer with finite thickness. While 
charge current is independent of position in ferromagnet, 
spin current is dependent on position and is not con- 
served. In fact, spin current at the left and right inter- 
faces have opposite sign and in the middle of ferromag- 
netic layer, the spin current vanishes. We have plotted 
spin current at the left interface between superconductor 
and ferromagnetic layer in Figs. 4[5 and Fig(8] At this 
point, there are three components of spin current. 

For charge current, we see that by increasing the thick- 
ness of ferromagnetic layer the current decreases with 
oscillation because of decreasing of quantum coherency, 
but for the spin current the situation is quite different. 
As a result of the exchange field, the thicker layer has a 
greater effect on spin of quasiparticles and thus a larger 
spin current. As seen in FigJ2J for a = 7r/2, the results 
are independent of the exchange field. Thus, charge cur- 
rent in this case is only dependent on thickness of the 
ferromagnetic layer and spin current is absent. Regard- 
ing spin currents, all components of spin current vanish 
at a = and a = ir/2. As shown in FigHJ spin current is 
strongly dependent on the exchange field, so for a weak 
ferromagnet (for example h = 0.5A), spin current is neg- 
ligibly small. For some values of the phase difference, 
spin current exists in the absence of charge current as 
seen from Figs. [4] and [5] Also, there is local maximum or 
minimum in spin current in <p = or <p — w. This can be 



explained as follows. Charge currents are generally odd 
function but spin currents are even function of the phase 
difference, because of the symmetry under the time re- 
versal. Therefore, the slope of the spin current becomes 
zero at <p = and = it. 

The critical currents as a function of h and I are 
shown in FigjS] and FigjT] for different misorientations 
a = 0, 7r/4 and tt/2. Also, critical charge and spin cur- 
rents as the function of a have been plotted in FigjH] 
Here, the critical charge and spin currents are defined 
as j ce = Max^j e (4>) = j e (4>*) and j s = j s {<j>*), respec- 
tively. The 0-7T transition as a function of h for large I 
and small a are shown in FigJ6] In RefpS it is shown 
that 0-7r transition cannot occur at a = tt/2 for a S func- 
tion type ferromagnetic barrier. We find that this is the 
case even in the presence of a ferromagnetic layer with 
finite thickness instead of a S function type barrier. This 
is also seen from the dependence of the Josephson cur- 
rent on I which shows the 0-n transition for large h and 
small a as seen in FigjTJ As shown in FigjHJ for finite 
thickness of ferromagnetic layer, 0-n transition can be 
found with changing the misorinetation angle. The 0-ir 
transition is accompanied by the jump of the spin, cur- 
rent due to the phase jump at the 0-7T transition.^ Even 
when charge current shows a monotonous dependence on 
a, spin current can show non-monotonous a dependence. 



V. CONCLUSIONS 

In this paper, we have investigated transport proper- 
ties in triplet SFS Josephson junction and clarified the 
effects of the thickness of ferromagnet, the exchange field, 
and misorientation angle between exchange field of ferro- 
magnet and (i-vector of superconductors on the Joseph- 
son charge and spin currents. We found that 0-n tran- 
sition can occur except for the case that the exchange 
field and d-vector are in nearly perpendicular configu- 
ration. Also, we showed how spin current flows due to 
misorientation between the exchange field and d-vector 
and found that it is absent when two vectors are parallel 
or perpendicular to each other. 
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